 function b = lagrange(a, x, y)
	% valoare polinom Lagrange in a
	%Intrari:
	% a = abscisa in care se cere polinomul
	% x = abscisele celor n+1 puncte
	% y = ordonatele celor n+1 puncte
	%Iesiri:valoare polinom interpolare in a
	n = length(x); % numar puncte
 	b = 0;
 	for i = 1 : n
   		produs = y(i);
   		for j = 1 : n
     			if i ~= j
        			produs = produs * ( a - x( j ) ) / ( x( i ) - x( j ) );
     			endif
   		endfor
   		b = b + produs;
 	endfor

endfunction
